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Abstract 

Recently a new model of galactic gravitational field, based on ordinary General Relativity has 
been proposed by Cooperstock and Tieu in which no exotic dark matter is needed to fit the 
observed rotation curve to a reasonable ordinary matter distribution. We argue that in this model 
' the gravitational field is generated not only by the galaxy matter, but by a thin, singular disk as 

o ' 

CN ' well. The model should therefore be considered unphysical. 
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I. INTRODUCTION 

In [lj the authors derived a set of non-linear equations which describes a gravitationally 
bound rotating cloud of dust in General Relativity. Basing on their result, the authors pro- 
pose a model of galaxy to which a flat rotation curve can be fitted without dark, unobservable 
matter. The result is then applied to several galaxies, including the Milky Way. 

We shall argue that 

1. the set of equations they derive is essentialy correct, though does not posses non- 
trivial, asymptotically flat solutions, 

2. the model of galaxy they propose is questionable because it posseses an additional 
source of gravitational field in the form of a rotating flat disk at z = 0. 



II. GRAVITATION ALY BOUND CLOUD OF DUST IN GR 

In GR a gravitationally bound cloud of dust is described by the stress-energy tensor 
T^ u = p u^u u such that is the dust four-velocity vector and p its density (see for example 
P]). The metric satisfies the Einstein equations 

G^[g af3 } = 8TrGpu»u» (1) 

where the Einstein tensor field is to be considered a functional of the metric field g a p (we 
set c = 1 throughout the paper). Note that the Biachi identities immediately imply that 
the flow of the dust is geodesic and source-free 

p/V/ = (2) 
V„(pu") = 0. (3) 

Let i] denote the flat Minkowski metric. The conventional way to simplify this highly 
non-linear set of equations is to perform a perturbative expantion of some of the quantities 
in G: 

9ap = Vap + G g§ + G 2 + 0(G 3 ) (4) 
u a = uf 0) + G ufo + G 2 uf 2) + 0(G 3 ) (5) 
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and obtain the following equation 

G ■ G'^[g m ] + G 2 ■ (G'^[g^} + ^G"^[g {1 \ + 0(G 3 ) = 
= 8nGp(u^ + Gvfc)(vfa + Gufa) + 0(G 3 ) 

(G'^ u [ • ] and G"^ u [ ■ , • ] denotes, respectively, the first and second functional derivative of 
G^ u [g a p\ for g a p = rjap). We can now attempt to solve these equations order by order. 

This approach, however, fails to reproduce the set of equations from This is due to 
the fact that the exact solution families approximated there depend on G via yG, i- e. their 
dependence on the coupling constant is non-analytic. Since we want to follow the approach 
of we shall expand (JTJ in \^G instead of G: 



and therefore 



gap = Vap + G 1 ' 2 g§ + G + 0(G 3 / 2 ) (6) 
u a = uf 0) + G l/2 uf 1} + G uf 2) + 0(G 3 / 2 ) (7) 



G i/2 . g^M] + G ( G'^lg®] + ^G'^lgW, g {1) ] ) + 0(G 3/2 ) 



= 8vrGp^ 0) ^ 0) + O(G 3 / 2 ). 
Rewritten order by order this is equivalent to 



G"H<?«] = (8) 
G'^[g {2) ] + l -G"^[g^\g^} = Snpu^u^. (9) 

The first equation is the linearized vacuum Einstein equation. If we impose additionally 
the Lorentz-De Donder gauge condition (see jsl), (JHJ) is equivalent to 

□<!) = (10) 

(□ being the flat metric rj d'Alembert operator). In particular, for stationary problems 

A^ = 0. (11) 

Note that for an asymptotically flat metric g^K — ► for r — > oo, which together with (fTTj) 
implies that g^ = 0. In other words, no asymptotically flat solutions are possible with the 
metric depending on yG in the lowest order. 



III. COOPERSTOCK AND TIEU GALACTIC ROTATION CURVE MODELS 



Cooperstock and Tieu assumed the metric to be of the form of 

g = - e u (dz 2 + dr 2 ) - r 2 d<p 2 + (dt - Nd<p) 2 (12) 

with 

N = r d - (13) 
$ = e - fc|z| Jo(kr), (14) 



both $ and N being of the order of yG. It is straightforward to verify that this metric 
fails to satisfy (JHJ), as Ag^ is not equal to zero, but has a distributional source proportional 
to S(z). (|T2*|) is not, therefore, a well-defined, global metric with rotating dust as the only 
source of gravity. 

In fact, it does not satisfy the first of the two equations derived by Cooperstock and Tieu: 

N 

N rr + N zz - = (15) 

r 

N 2 + N 2 

= 8nG P (16) 

As before, due to the presence of \z\ in (fT4*j) the first equation has a distributional singularity 
at z = rather than on the right-hand side. In the following section we shall prove that 
this results in presence of an additional, distributional term in the stress-energy tensor T^ v . 



IV. KOMAR INTEGRAL 

The shortest way to prove that the galaxy model in Q has a sheet of matter at z = is 
to consider the Komar integral of the timelike Killing vector Define r = g(-jL, •)• It is 
easy to prove in general that if V^Tj, + V v t^ = 0, the following identity holds: 

d * dr = - r a e^pv dx u A dx* 3 A dx a (17) 
3 

(e is the volume 4-form; if the metric is vacuum the identity reduces to the Komar current 
conservation formula). We can integrate it over an arbitrary 3-volume V: 

f *dr = f d * dr = ^— f (2T a »T a - Tr^e^ dx v A dx" A dx a . (18) 

JdV JV 3 Jy 
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Take V to be the cylinder defined by the following conditions: r < R, —a < z < a, 
t = const, for given positive constants a and R. The boundary of this volume consists of a 
side 2-surface r = R and two circle-shaped 2-surfaces z = a and z = —a. Let us take the 
limit of a — > 0, i. e. the flat cylinder limit. If T^ v had no singularities but consisted solely 
of rotating dust, this limit would be zero, since the volume of the cylinder goes to zero as 
we shrink it. It is, however, not the case in the metric (fT2"j) . Indeed, the side surface integral 
of (fTHj) does converge to 0, but the top and bottom side integrals I t and I b do not: 
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I t + I b = dy dr 



N dN 



2t , f R N dN 
dip / dr 



r dz 



(19) 



r dz 

With iV given by (j!3|) and (fTlj) I t and 1^ are in fact equal and therefore 

r R N BN f R 
I t + I b = A7r dr — — -> -Ank 3 / dr r J' Q {krf ^ 0. (20) 

JO r & Z z=a JO 

The reason is, once again, that ^ has a discontinuity at z = due to the presence of \z\ in 
its definition. 

Note that the same argument can be applied to the axial Killing vector as well. This 
indicates that the disk has both mass and angular momentum, simmilarily to the singularity 
in the exact solution of van Stockum Q]. 

The solution should therefore be interpreted as possesing an additional, distributional 
source of gravity at z = and the resulting discontinuity of the metric's first derivative as 
a physical one rather than a coordinate system artifact, as is claimed in 
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